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Magnetic helicity is nearly conserved and its evolution equation provides a 
dynamical feedback on the alpha effect that is distinct from the conventional 
algebraic alpha quenching. The seriousness of this dynamical alpha quenching 
is particularly evident in the case of closed or periodic boxes. The explicit 
connection with catastrophic alpha quenching is reviewed and the alleviating 
effects of magnetic and current helicity fluxes are discussed. 



1 Introduction 

Let us begin by defining dynamos and helicity. Dynamos are a class of ve- 
locity fields that allow a weak seed magnetic field to be amplified until some 
saturation process sets in. Mathematically, this is described by exponentially 
growing solutions of the induction equation. Simulations have shown that 
any sufficiently complex flow field can act as a dynamo if the resistivity is 
below a certain threshold. It is in principle not even necessary that the flow 
is three-dimensional, only the magnetic field must be three-dimensional be- 
cause otherwise one of several antidynamo theorems apply (Cowling 1934, 
Zeldovich 1957). 

Helicity, on the other hand, quantifies the swirl in a vector field. There 
is kinetic helicity, which describes the degree to which vortex lines follow a 
screw-like pattern, and it is positive for right-handed screws. Examples of 
helical flows are the highs and lows on the weather map. For both highs 
and lows the kinetic helicity has the same sign and is negative (positive) in 
the northern (southern) hemisphere. For example, in an atmospheric low, air 
flows inward, i.e. toward the core of the vortex, and down to the bottom of 
the atmosphere, but the Coriolis force makes it spin anti-clockwise, causing 
left-handed spiraling motions and hence negative helicity. 

A connection between helicity and dynamos has been established already 
quite some time ago when Steenbeck et al. (1966) calculated the now famous a 
effect in mean field dynamo theory and explained its connection with kinetic 
helicity. In this paper we are not so much concerned with kinetic helicity, 
but mostly with the magnetic and current helicities. Quantifying the swirl 
of magnetic field lines has diagnostic significance, because magnetic helicity 
is a topological invariant of the ideal (non-resistive) equations. Especially 
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in the solar community the diagnostic properties of magnetic helicity have 
been exploited extensively over the past decade. However, the use of magnetic 
helicity as a prognostic quantity for understanding the governing nonlinearity 
of a effect dynamos has only recently been noted in connection with the 
magnetic helicity constraint (Brandenburg 2001, hereafter referred to as B01). 

We should emphasize from the beginning that dynamos do not have to 
have helicity. The small scale dynamo of Kazantsev (1968) is an example of 
a dynamo that works even without helicity. Nonhelical dynamos are gener- 
ally harder to excite than helical dynamos, but both can generate fields of 
appreciable strength if the magnetic Reynolds number is large. The stretch- 
twist-fold dynamo also operates with twist (as the name suggests!), but the 
orientation of twist can be random, so the net helicity can be zero. Simula- 
tions have shown that even with zero helicity density everywhere, dynamos 
can work (Hughes et al. 1996). 

It is also possible to generate magnetic fields of large scale once there 
is strong shear, even if there is no helicity (Vishniac & Brandenburg 1997). 
This case is very much a topic of current research. One of the possibilities 
is is the so-called shear-current effect (Rogachevskii & Kleeorin 2003, 2004), 
but such dynamos still produce helical large scale magnetic fields. There 
is also the possibility of an intrinsically nonlinear dynamo operating with 
magnetic helicity flux alone (Vishniac & Cho 2001). Thus, it is not necessarily 
clear that large scale dynamos have to work with kinetic helicity and the 
corresponding a effect. However, there is as yet no convincing example of a 
dynamo without the involvement of kinetic helicity that generates large scale 
magnetic fields with a degree of coherence that is similar to that observed in 
stars and in galaxies, e.g. cyclically migrating magnetic fields in the sun and 
grand magnetic spirals in some nearby galaxies. Such fields can potentially 
be generated by dynamos with an a effect, as has been shown in many papers 
over the past 40 years; see Chapters 2, 4, and 10. 

There is however a major problem with a effect dynamos; sec Branden- 
burg (2003), Brandenburg & Subramanian (2004) for recent reviews on the 
issue. The degree of severity depends on the nature of the problem. It is 
most severe in the case of a homogeneous a effect in a periodic box, which 
is also when the problem shows up most pronouncedly. Cattaneo & Hughes 
(1996) found that the a effect is quenched to resistively small values once the 
mean field becomes a fraction of the cquipartition field strength. In response 
to such difficulties three different approaches have been pursued. The most 
practical one is to simply ignore the problem and the proceed as if we can 
still use the a effect with a quenching that only sets in at equipartition field 
strengths. This can partially be justified by the apparent success in applying 
this theory; see the recent reviews by Beck et al. (1996), Kulsrud (1999), and 
Widrow (2002). The second approach is to resort to direct three dimensional 
simulations of the turbulence in such astrophysical bodies. In the solar com- 
munity this approach has been pioneered by Gilman (1983) and Glatzmaier 
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(1985), and more recently by Brun et al. (2004). The third approach is a 
combination of the first two, i.e. to use direct simulations of problems where 
mean field theory should give a definitive answer. This is also the approach 
taken in the present work. The hope is ultimately to find guidance toward 
a revised mean field theory and to test it quantitatively. A lot of progress 
has already been made which led to the suggestion that only a dynamical 
(i.e. explicitly time dependent) theory of a quenching is compatible with the 
simulation results. In the present paper we review some of the simulations 
that have led to this revised understanding of mean field theory. 

The dynamical quenching theory is now quite successful in reproducing 
the results from simulations in a closed or periodic box with and without 
shear. In these cases super-equipartition fields are possible, but only after 
a resistively long time scale. In the case of an open box without shear the 
dynamical quenching theory is also successful in reproducing the results of 
simulations, but here the root mean field strength decreases with increas- 
ing magnetic Reynolds number, suggesting that such a dynamo is unimpor- 
tant for astrophysical applications. Open boxes with shear appear now quite 
promising, but the theory is still incomplete and, not surprisingly, there are 
discrepancies between simulations and theory. In fact, it is quite possible 
that it is not even the a effect that is important for large scale field regenera- 
tion. Alternatives include the shear-current effect of Rogachevskii & Kleeorin 
(2003, 2004) and the Vishniac & Cho (2001) magnetic and current helicity 
fluxes. In both cases strong helicity fluxes are predicted by the theory and 
such fluxes are certainly also confirmed observationally for the sun (Berger 
& Ruzmaikin 2000, DeVore 2000, Chae 2000, Low 2001). For the galaxy the 
issue of magnetic helicity is still very much in its infancy, but some first 
attempts in this direction are already being discussed (Shukurov 2004). 

2 Dynamos in a periodic box 

To avoid the impression that all dynamos have to have helicity, we begin by 
commenting explicitly on dynamos that do not have net kinetic helicity, i.e. 
\{u ■ u)\/(k[(u 2 )) <C 1, where fcf is the wavenumber of the forcing (corre- 
sponding to the energy carrying scale). Unless the flow also possesses some 
large scale shear flow (discussed separately in Sect. 14 .51 below ). such dynamos 
are referred to as small scale dynamos. The statement made in the introduc- 
tion that any sufficiently complex flow field can act as a dynamo is really 
only based on experience, and the statement may need to be qualified for 
small scale dynamos. Indeed, whether or not turbulent small scale dynamos 
work in stars where the magnetic Prandtl numbers are small (PrM w 10 -4 ) 
is unclear (Schekochihin et al. 2004, Boldyrev & Cattaneo 2004). Simulations 
suggest that the critical magnetic Reynolds numbers increase with decreasing 
magnetic Prandtl number like i? m ,crit ~ 35Pi'm (Haugcn ct al. 2004). 
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Throughout the rest of this review, we want to focus attention on large 
scale dynamos. This is where magnetic helicity plays an important role. Be- 
fore we explain why in a periodic box nonlinear dynamos operate only on 
a resistively slow time scale, it may be useful to illustrate the problem with 
some numerical facts. 

In the simulations of B01 the flow was forced at an intermediate wavenum- 
ber, k ~ kf — 5, while the smallest wavenumber in the computational domain 
corresponds to A; = fci = 1. The kinetic energy spectrum peaks at k « ki, 
which is therefore also the wavenumber of the energy carrying scale. The 
turbulence is nearly fully helical with (u> ■ u) / (kf(u 2 )) 0.7. ..0.9. The initial 
field is random, but as time goes on it develops a large scale component at 
wavenumber k ps k\ = 1; see Fig. ^ In Fig. [3 we show the evolution of the 
magnetic energy of the mean field from the same simulation. 1 Here the mean 
field is defined as two-dimensional averages over planes perpendicular to the 
direction in which the mean field varies. There are of course three such direc- 
tions, but there is usually only one direction for which there is a significant 
mean field. 

While the saturation field strength increases with increasing magnetic 
Reynolds number, the time scale on which this nonlinear dynamo saturates 
increases too. To avoid misunderstandings, it is important to emphasize that 
this result applies only when we are in the nonlinear regime and when the 
flows are helical. 

In turbulence one is used to situations where the microscopic values of 
viscosity v and magnetic diffusivity r\ do not matter in the sense that, for 
almost all practical purposes, they are superseded by turbulent effective val- 
ues, ft and r/t, respectively. This is because in turbulence there is spectral 
energy all the way down to the viscous/resistive length scale, (t^t) 1 / 2 , where 
r is the turnover time. 2 Thus, even when v is very small, the rate of vis- 
cous dissipation, {2vpS 2 ), is in general finite (S is the trace-less rate of strain 
tensor). Likewise, even when r\ is very small, the rate of Joule dissipation, 
r]fj, (J 2 ), is in general finite (^o is the magnetic permeability). This is because 
the current density diverges with decreasing n like \ J\ ~ rj -1 ^ 2 , so the energy 
dissipation stays finite and asymptoticly independent of how small r\ is. The 
trouble is that the value of magnetic helicity dissipation is proportional to 
?y( J ■ B) (see below), and in the limit rj — > we have tj( J ■ B) — > rj 1 ^ 2 — > 0, so 
resistive magnetic helicity dissipation becomes impossible in the limit of large 
R m . In the following section we derive and discuss the evolution equation for 
magnetic helicity. 

1 Here the time unit is [t] = (c s fci) _1 , where c s is the isothermal speed of sound, 
and the magnetic field is measured in units of [B] = */]Iopoc B . 

2 The turnover time at the wavenumber k is («^fc) _1 . Using Kolmogorov scal- 
ing, Uk ~ fe -1 / 3 , one finds the familiar formula k v — kiB^l A , where kf is the 
wavenumber of the energy carrying eddies. 
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Fig. 1. Cross-sections of B x (x,y,0) for Run 3 of B01 at different times showing 
the gradual build-up of the large scale magnetic field after t = 300. The diffusive 
time scale for this run is (r/kl) -1 = 500. Dark (light) corresponds to negative 
(positive) values. Each image is scaled with respect to its min and max values. 
The final state corresponds to the second eigenfunction given in Eq. (ESI, but 
with some smaller scale turbulence superimposed. [Adapted from Brandenburg & 
Subramanian (2004).] 



3 Magnetic helicity evolution 

3.1 The two-scale property of helical turbulence 

Usually in mean field dynamo theory one talks about the two-scale assump- 
tion made in order to derive the mean field equations (e.g. Moffatt 1978, 
Krause & Radler 1980). This has to do with the fact that higher order deriva- 
tives in the mean field equation can only be neglected when the mean field is 
sufficiently smooth. Here, instead, we use the two-scale properties of helical 
turbulence as demonstrated in the previous section. These properties emerge 
automatically when the size of the whole body is at least several times larger 
than the scale of the turbulent eddies. As Fig.^shows explicitly, a large scale 
field (wavenumber k\) emerges in addition to the forcing scale (wavenumber 
k{ > fci). 

In this section we discuss the magnetic helicity equation and use it to- 
gether with the two-scale property of helical turbulence to derive the so-called 
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Fig. 2. Evolution of {b' 2 ) for Runs 1-3 and 5, 6 (dashed lines). The magnetic 
Reynolds numbers are _R m = u TraB /rjkf = 2.4, 6, 18, 100, and 16, respectively; 
see B01. The solid lines denote the solution of the associated mean field dynamo 
problem where both a and and turbulent diffusivity r/t are quenched in a magnetic 
Reynolds number dependent fashion. [Adapted from B01.] 

magnetic helicity constraint that allows the result of Fig. [21 to be understood 
quantitatively. 

3.2 Definition of helicity 

The helicity of any solenoidal vector field /, i.e. with V • / = 0, is defined 
as the volume integral of / dotted with its inverse curl, i.e. curl / = g. As 
pointed out by Moffatt (1969), the helicity quantifies the topological linkage 
between tubes in which / is non- vanishing. In the following the linkage aspect 
of helicity will not be utilized, but rather the mathematical evolution equation 
that the helicity obeys (see the next section). However, the calculation of 
g is problematic because it involves a gauge ambiguity in that the curl of 
g' = g + also gives the same / = curlg'. 

In the special case of periodic boundary conditions or for h ■ f = on 
the boundaries, where h is the normal vector on the boundary, the helicity 
is actually gauge- invariant, because 
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J f g'dV = J f gdV + J f-Vc^dV 

= [ fgdV- [ <f>VfdV= [ fgdV, (1) 

where we have used V • /. Since the magnetic field is divergence free, the 
magnetic helicity, J B ■ cm\^ 1 BdV is gauge invariant. For other boundary 
conditions this is unfortunately not the case. 

For vector fields whose inverse curl is a physically meaningful quantity, 
such as the vorticity u>, whose inverse curl is the velocity u, the gauge ques- 
tion never arises. In this and similar cases the helicity density, u) ■ u in 
this case, is physically meaningful. Other examples are the cross helicity, 
J B ■ curl _1 o;dF, which describes the linkage between magnetic flux tubes 
and vortex tubes, and the current helicity, J J ■ curl 1 J dV, which quantifies 
the linkage of current tubes. In these two cases it is natural to use curl _1 = B 
and curl = u. For the magnetic field one can define the magnetic vector 
potential, curl _1 S = A, but A is not a physically meaningful quantity and 
hence the magnetic helicity, 

H = [ A-BdV= (A-B)V (2) 
Jv 

is gauge-dependent, unless the boundaries of the volume V are periodic or 
perfectly conducting. Here and below, angular brackets denote volume av- 
erages. Occasionally, however, we simply refer to {A ■ B) as the magnetic 
helicity, but this is strictly speaking only the magnetic helicity per unit vol- 
ume. 

In the following section we derive the evolution equation for (A ■ B) and 
focus first on the case where the boundary conditions are indeed periodic, so 
(A ■ B) is gauge- invariant. 

3.3 Derivation of the magnetic helicity equation 

The homogeneous Maxwell equations are 

— = -VxE, V-B = 0. (3) 

Expressing this in terms of the magnetic vector potential, A, where B — 
V x A, we have 

8 A 

M=- E -^ (4) 
where cf> is the scalar potential. Dotting Eqs © and with A and B, 
respectively, and adding them, we have 

^-(A ■ B) = —2E ■ B -V ■ (E x A + <j>B). (5) 
ot 
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Here, A ■ B is the magnetic helicity density, but since it is not gauge invari- 
ant (see below) it is not a physically meaningful quantity. After integrating 
Eq. (JSJ over a periodic volume, the divergence term does not contribute. Fur- 
thermore, using Ohm's law, E = —U x B + TjfioJ, where J = V x B//j,q is 
the current density, we have 

^ t (A-B) = -2 Wo (J-B), (6) 

i.e. the magnetic helicity, (A ■ B), changes only at a rate that is proportional 
to r)(J -B). (Here and elsewhere, angular brackets denote volume averaging.) 
As discussed in the previous section, this rate converges to zero in the large 
i? m limit. Here, angular brackets denote volume averages, i.e. (A ■ B) = 
±J v A-BdV. 

We recall that for periodic boundary conditions, (A-B) is invariant under 
the transformation A — > A' = A + V/l, which does not change the value of 
B' = V x A' = V x A = B. Here, A is a gauge potential. Thus, for periodic 
boundary conditions, (A-B) is a physically meaningful quantity. The same is 
also true for perfectly conducting boundaries (see Brandenburg & Dobler 2002 
for corresponding simulations). For open boundaries, however, J v A- B dV is 
not gauge invariant, but one can derive a gauge-invariant relative magnetic 
helicity (Berger & Field 1984). 



3.4 The magnetic helicity constraint 

A very simple argument can be made to explain the saturation level and the 
resistively slow saturation behavior observed in Fig. The only assumption 
is that the turbulence is helical, i.e. (u) ■ u) ^ 0, where us is the vorticity, and 
that this introduces current helicity (j -b), at the same scale and of the same 
sign as the kinetic helicity. Here we have split the field into large and small 
scale fields, i.e. B = B + b and hence also J = J + j and A = A + a. 

The first remarkable thing to note is that, even though we are dealing 
with helical dynamos, there is no net current helicity in the steady state, i.e. 

(J-B)= 0; (7) 

see Eq. However, using the decomposition into large and small scale 
fields, we can write 

(J-B) = (J-B) + (j-b)=0, (8) 

so we have 

(J-B) = -(j-b) (9) 
in the steady state. We now introduce the approximations 3 

3 Here and elsewhere we use units where fio — 1 or, following R. Blandford (private 
communication), we use units in which pi is one quarter. 
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(J-B)^k m (B 2 ) and (j ■ b) w ~k f (b 2 ) (10) 

where k m and kf are the typical wavenumbers of the mean and fluctuating 
fields, respectively. These approximations are only valid for fully helical tur- 
bulence, but can easily be generalized to the case of fractional helicity; see 
Sect. O and Blackman & Brandenburg (2002, hereafter BB02). We have 
furthermore assumed that the sign of the kinetic helicity is negative, as is 
the case in the northern hemisphere of the sun, for example. (The case of 
positive kinetic helicity is straightforward; see below.) The wavenumber k{ of 
the fluctuating field is for all practical purposes equal to the wavenumber of 
the forcing function. (In more general situations, such as convection or shear 
flow turbulence, k{ would be the wavenumber of the energy carrying eddies.) 
We also note that for large values of the magnetic Reynolds number, J? m , the 
kf factor in Eq. ljTU|l gets attenuated by an rU 4 factor (BB02). On the other 
hand, the wavenumber of the mean field is in practice the wavenumber of the 
box, i.e. k m = k\. Inserting now Eq. (fTUfl into yields 

(S 2 > = ^(b 2 ), (11) 

i.e. the energy of the mean field can exceed the energy of the fluctuating 
field - in contrast to earlier expectations (e.g. Vainshtein & Cattaneo 1992, 
Kulsrud & Andersen 1992, Gruzinov & Diamond 1994). Indeed, in the two- 
dimensional case there is an exact result due to Zeldovich (1957), 

(B 2 ) = i? m 1 (6 2 ) (2-dimensional case). (12) 

This result has also be derived in three dimensions using first order smoothing 
(Krause & Radler 1980), but it is important to realize that this result can 
break down in the nonlinear case in three dimensions, where Eq. (|llfl is in 
good agreement with the simulations results. However, the assumption of 
periodic or closed boundaries is an essential one. We return to the more 
general case in Sects l4~4ll4. 51 

The time dependence near the saturated state can be approximated by 
using 

( J • B) « k 2 m (A ■ B) and ( j ■ b) w k\ (a-b). (13) 

These equations are still valid in the case of fractional helicity (BB02). Only 
the two-scale assumption is required. Near saturation, 

\(A-B)\= (£)V b >< ( 14 ) 

i.e. \(A ■ B)\ ^S> | (a ■ b)\ and so we can neglect (a-b), and the magnetic 
helicity equation @ becomes therefore an approximate evolution equation 
for the magnetic helicity of the mean field, 
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dV 

or, by using Eq. i fTUj l 



d 'A-B) = -2 Wo (J ■ B) - 2 Wa (j ■ b), (15) 



-i d /-n 2 \ _ o-i. m\ , ..2.. 



fc m i -(S > = -2nk m (B ) + 2j7fe f (6^). (16) 
Note the plus sign in front of the (b 2 ) term resulting from Eq. I|10[) . The plus 

2 

sign leads to growth while the minus sign in front of the (B ) term leads to 
saturation (but both terms are proportional to the microscopic value of ij; 
see below) . Once the small scale field has saturated, which will happen after 
a few dynamical time scales such that (b 2 ) « B 2 q = fio(pu 2 }, the large scale 
field will continue to evolve slowly according to 



(B 2 ) = ^-(b\ 

u2 



-2»jfe4(t-t s »0 



(17) 



where £ sat is the time at which (b ) has reached approximate saturation. 
In practice, t sa ± can be determined such that Eq. I|17|l describes the simula- 
tion data best. We refer to Eq. (|17fl as the magnetic helicity constraint. The 
agreement between this and the actual simulations (Fig. |3J) is quite striking. 

The significance of this remarkable and simple equation and the almost 
perfect agreement with simulations is that the constraint can be extrapolated 
to large values of R m where it would provide a benchmark, against which all 
analytic dynamo theories, when subjected to the same periodic boundary 
conditions, should be compared to. In particular the late saturation behavior 
should be equally slow. We return to this in Sect. [51 

An important question is whether anything can be learned about stars 
and galaxies. Before this can be addressed, we need to understand the differ- 
ences between dynamos in real astrophysical bodies and dynamos in periodic 
domains. 



4 What do stars and galaxies do differently? 

We begin with a discussion of fractional helicity, shear and other effects that 
cause the magnetic helicity to be reduced. We then address the possibility of 
helicity fluxes through boundaries, which can alleviate the helicity constraint 
(Blackman & Field 2000). 

4.1 Fractional helicity 

When the turbulence is no longer fully helical, Eq. (|10|) is no longer valid and 
needs to be generalized to 

(J • B) = e nl k m (B 2 ) and (j ■ b) - -e f fc £ (b 2 ), (18) 
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Fig. 3. Late saturation phase of fully helical turbulent dynamos for three differ- 
ent values of the magnetic Reynolds number: i? m = u ITas /rik{ = 2.4, 6, and 18 for 
Runs 1, 2, and 3 respectively; see B01. The mean magnetic field, B, is normalized 
with respect to the equipartition value, B cq — y/iopoMrms, and time is normalized 
with respect to the kinematic growth rate, A. The dotted lines represent the fit 
formula l|17|l which tracks the simulation results rather well. [Adapted from Bran- 
denburg et al. (2003).] 

where e m < 1 and ef < 1 are coefficients denoting the degree to which the 
mean and fluctuating fields are helical. Equation (jl3|l is still approximately 
valid in the fractionally helical case. 

Maron & Blackman (2002) found that there is a certain threshold of £f 
below which the large scale dynamo effect stops working. Qualitatively, this 
could be understood by noting that the large scale magnetic field comes from 
the helical part of the flow, so the velocity field can be though of as having 
a helical and a nonhclical component, i.e. 

U = U hcl + U nohcl . (19) 

However, the dynamo effect has to compete with turbulent diffusion which 
comes from both the helical and the nonhelical parts of the flow. Thus, when 
|fnohei| becomes too large compared with \Uhd\ the large scale dynamo effect 
will stop working. 
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Although we have not yet discussed mean field theory we may note that 
the value of the threshold can be understood quantitatively (Brandenburg et 
al. 2002, hereafter BDS02) and one finds that large scale dynamo action is 
only possible when 

k 

ef > (for large scale dynamos). (20) 

fc f 

In many three-dimensional turbulence simulations or in astrophysical bodies, 
this threshold criterion may not be satisfied, and so mean field dynamo of 
the type described above (a 2 dynamo) may not be excited. If there is shear, 
however, this criterion will be modified to 

et > e m ^, (21) 

where ef is the degree to which the large scale field is helical. In dynamos 
with strong shear, | ef | may be very small, making mean field dynamo action 
in fractionally helical flows more likely. This will be discussed in the next 
section. 



4.2 Dynamos with shear 

In the presence of shear, the streamwise component of the field can be am- 
plified by winding up the poloidal (cross-stream) field. Again, the resulting 
saturation field strength can be estimated based on magnetic helicity conser- 
vation arguments. 

Note first of all that for closed or periodic domains, Eq. (JSJ) is still valid 
and therefore (J • B) = —(j ■ b) in the steady state. 4 However, while (j ■ b) 
still depends on the helicity of the small scale field, the corresponding value 
of ( J ■ B) no longer provides such a stringent bound on (B ) as before. This 
is because shear can amplify the toroidal field independently of any magnetic 

2 

helicity considerations. The component of B that is amplified by shear is 
nonhelical and so we have 

e m = |(7-B)|/(fc m (B 2 )) «1 (22) 

(or at least |e m | <C 1 when the helicity of the forcing is negative and e m 
therefore negative). The value of e m is proportional to the ratio of poloidal 
to toroidal field 

e m ^±2(B 2 pol )^/(Biy/\ (23) 

4 This is because in the E ■ B term in the magnetic helicity equation the induction 
term, U x B, drops out after dotting with B. (For this reason, also ambipolar 
diffusion and the Hall effect do not change magnetic helicity conservation.) 
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Fig. 4. Growth of poloidal and toroidal magnetic fields on a logarithmic scale 
(upper panel), and product of poloidal and toroidal magnetic fields on a linear 
scale. For the fit we have used k\ — 2, B cq — 0.035, and eo = 1-3. [Adapted from 
BBS01.] 



where the numerical pre-factor can be different for different examples. 5 With 
these preparations the magnetic helicity constraint can be generalized to 



2B™?B? mB ~ fcf "~ 2 " 



'pol -"tor 



1 - e 



-2»)/4(i-t aat ) 



(24) 



This form of the constraint was proposed and confirmed using three-dimensional 
simulations of forced helical turbulence with large scale shear (Brandenburg 
et al. 2001, hereafter BBS01); see also Fig. EJ 

The main conclusion to be drawn from this is that the magnetic helicity 
constraint is still valid in the presence of shear, i.e. the timescale of saturation 

5 Take as an example B(z) — (B po \, B t0 r, 0) T = (e cos k\z, sin kiz,0) T for e«l, 
so {B ) ~ 1/2 and {J ■ B) ~ eki and therefore e m = 2e. 
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is still controlled by the microscopic magnetic diffusivity. The only difference 
is that stronger field strengths are now possible. 

Another interesting aspect is that dynamos with shear allow for oscilla- 
tory solutions of the magnetic field. This is expected from mean field dynamo 
theory (Steenbeck & Krause 1969a,b), but it is also borne out by simulations 
(BBS01). The main result is that the resulting cycle frequency seems to scale 
with the microscopic magnetic diffusivity, not the turbulent magnetic dif- 
fusivity. This confirms again that in a closed domain the magnetic helicity 
constraint plays a crucial role in controlling the timescale of nonlinear dy- 
namos. 

4.3 Hall effect dynamos 

In recent years the importance of the Hall effect has been emphasized by a 
number of groups, especially in applications to protostellar accretion discs 
(Balbus & Terquem 2001). The hall effect can lead to strong nonlinear steep- 
ening of field gradients (Vainshtein et al. 2000), and is therefore important 
for fast reconnection (e.g. Rogers et al. 2001), which in turn is relevant for 
neutron stars (Hollcrbach & Riidiger 2004). Nevertheless, since magnetic he- 
licity generation (and removal) is proportional to the dot product of electric 
and magnetic fields, and since the Hall current is proportional to J x B, the 
Hall term does not affect magnetic helicity conservation. Therefore the resis- 
tively limited saturation behavior should not be affected by the Hall term. 
Nevertheless, some degree of extra field amplification of the large scale field 
has been reported (Mininni et al. 2003), and it will be interesting to identify 
exactly the processes that led to this amplification. 

4.4 Magnetic helicity exchange across the equator or with depth 

The presence of an equator provides a source of magnetic helicity exchange 
between domains of negative helicity in the northern hemisphere (upper disc 
plane in an accretion disc) and positive helicity in the southern hemisphere 
(lower disc plane) . A similar situation can also arise in convection zones where 
the helicity is expected to change with depth (Yoshimura 1975). 

So far, simulations have not yet shown that the losses of small scale mag- 
netic fields are actually stronger than those of large scales fields. In Fig. [3] we 
show the saturation behavior of a system that is periodic, but the helicity of 
the forcing is modulated in the z-direction such that the sign of the kinetic 
helicity changes in the middle. One can therefore view this system as two 
subsystems with a boundary in between. This boundary would correspond to 
the equator in a star or the midplane in a disc. It can also model the change 
of sign of helicity at some depth in a convection zone. 

As far as the magnetic helicity constraint is concerned, the divergence 
term of current helicity flux is likely to be important when there is a bound- 
ary between two domains with different helicities. Naively, one might expect 
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Fig. 5. Evolution of the magnetic energy for a run with homogeneous forcing func- 
tion (solid line) and a forcing function whose helicity varies sinusoidally throughout 
the domain (dotted line) simulating the effects of equators at the two nodes of the 
sinusoidal helicity profile. [Adapted from Brandenburg et al. (2003).] 

there to be current helicity fluxes that are proportional to the current helic- 
ity gradient, analogous to Fick's diffusion law. These current helicity fluxes 
should be treated separately for large and small scale components of the field, 
so we introduce approximations to the current helicity fluxes from the mean 
and fluctuating fields as 

T m w -ThnVCU, « -iJfVQ. (25) 

The rate of magnetic helicity loss is here proportional to some turbulent 
diffusivity coefficient, 7j m or rji for the losses from mean or fluctuating parts, 
respectively. We assume that the small and large scale fields are maximally 
helical (or have known helicity fractions e m and ef) and have opposite signs 
of magnetic helicity at small and large scales. The details can be found in 
BDS02 and Blackman & Brandenburg (2003). The strength of this approach 
is that it is quite independent of mean field theory. 

We proceed analogously to the derivation of Eq. l(T7|) where we used the 
magnetic helicity equation © for a closed domain to estimate the time deriva- 
tive of the magnetic helicity of the mean field by neglecting the time derivative 
of the fluctuating field. This is a good approximation after the fluctuating 
field has reached saturation, i.e. t > i S at- Thus, we have 

k- 1 -^B 2 = -2 Vnl k m B 2 + 2 m k f b 2 , (26) 

where r) m — rji — r\ corresponds to the case of a closed domain. Note also 
that we have here ignored the volume integration, so we are dealing with 
horizontal averages that depend still on height and on time. 
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After the time when the small scale magnetic field saturates, i.e. when 
t > i S at7 we have (b 2 ) w constant. After that time, Eq. I|26|) can be solved to 
give 



for t > t sat . (27) 



(B 2 ) = (b 2 )^L [l _ e -2r )m fci( t -t sat) " 

This equation demonstrates three remarkable properties (Brandenburg et al. 
2003, Brandenburg & Subramanian 2004): 

Large scale helicity losses are needed (ry m > 77) to shorten the typical time 
scale. This is required to prevent resistively long cycle periods. 
However, the saturation amplitude is proportional to r){/r) m , so the large 
scale field becomes weaker as ?y m is increased. Thus, 
also small scale losses are needed to prevent the saturation amplitude 
from becoming too small. 

Future work can hopefully verify that these conditions are indeed obeyed 
by a working large scale dynamo. Simulations without shear have been unsuc- 
cessful to demonstrate that small scale losses are important (Brandenburg & 
Dobler 2001), but new simulations with shear now begin to show significant 
small scale losses of current helicity, an enhanced a effect (Brandenburg & 
Sandin 2004), and strong large scale dynamo action (see below). 



4.5 Open surfaces and shear 

The presence of an outer surface is in many respects similar to the presence of 
an equator. In both cases one expects magnetic and current helicity fluxes via 
the divergence term. A particularly instructive system is helical turbulence 
in an infinitely extended horizontal slab with stress-free boundary conditions 
and a vertical field condition, i.e. 

U X ,Z = Uy iZ = U Z = B X = By = 0. (28) 

Such simulations have been performed by Brandenburg & Dobler (2001) who 
found that a mean magnetic field is generated, similar to the case with pe- 
riodic boundary conditions, but that the energy of the mean magnetic field, 
2 

(B ), decreases with magnetic Reynolds number. Nevertheless, the energy of 
the total magnetic field, (B 2 ), does not decrease with increasing magnetic 

2 —1/2 

Reynolds number. Although they found that (B ) decreases only like R m , 
new simulations confirm that a proper scaling regime has not yet been reached 
and that the current data may well be compatible with an i?" 1 dependence; 
see Fig. EJ 

Clearly, an asymptotic decrease of the mean magnetic field must mean 
that the small scale dynamo does not work with such boundary conditions. 
Thus, the anticipated advantages of open boundary conditions are not borne 
out by this type of simulations. 
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Fig. 6. Dependence of the energy of the mean magnetic field on the magnetic 
Reynolds number for a run with open boundary conditions and no shear. 




Fig. 7. Left: A sketch of the solar angular velocity at low latitudes with spoke-like 
contours in the bulk of the convection zone merging gradually into uniform rotation 
in the radiative interior. The low latitude region, modeled in this paper, is indicated 
by thick lines. Right: Differential rotation in our cartesian model, with the equator 
being at the bottom, the surface to the right, the bottom of the convection zone 
to the left and mid-latitudes at the top. [Adapted from Brandenburg & Sandin 
(2004).] 

At this point we can mention some new simulations in a cartesian domain 
where differential rotation has been modeled as a region of the convection 
zone without explicitly allowing for convection; see Fig- El Instead, an exter- 
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Fig. 8. Visualization of the toroidal magnetic field during three different times 
during the growth and saturation for the run without kinetic helicity. 

rial forcing term has been applied that also drives the differential rotation. 
(Studies of the a effect have already been published; see Sect. 15.61 for details 
of the simulations and Sect. [S] for a discussion of the direct correspondence 
between the helicity constraint and the so-called catastrophic a quenching.) 
Here we briefly report on recent explicit dynamo simulations that have been 
carried out in this geometry. 

The size of the computational domain is ^tt x 2n x ^tt and the numerical 
resolution is 128 x 512 x 128 meshpoints. The magnetic Reynolds number 
based on the forcing wavenumber and the turbulent flow is around 80 and 
shear flow velocity exceeds the rms turbulent velocity by a factor of about 
5. We have carried out experiments with no helicity in the forcing (labeled 
by a — 0), as well as positive and negative helicity in the forcing (labeled by 
a < and a > 0, respectively); see Fig.|Hlfor a visualization of the run without 
kinetic helicity. We emphasize that no explicit a effect has been invoked. The 
labeling just reflects the fact that, in isotropic turbulence, negative kinetic 
helicity (as in the northern hemisphere of a star or the upper disc plane in 
galaxies) leads to a positive a effect, and vice versa. 

We characterize the relative strength of the mean field by the ratio q = 
(B )/{B 2 ), where overbars denote an average in the toroidal (y) direction; 
see Fig. El There are two surprising results emerging from this work. First, in 
the presence of shear rather strong mean fields can be generated, where up to 
70% of the energy can be in the mean field; see Fig. Second, even without 
any kinetic helicity in the flow there is strong large scale field generation. 
Obviously, this cannot be an aft dynamo in the usual sense. One possibility 
is the 6 x J effect, which emerged originally in the presence of the Coriolis 
force; see Radler (1969) and Krause & Radler (1980). In the present case 
with no Coriolis force, however, a S x J effect is possible even in the presence 
of shear alone, because the vorticity associated with the shear contributes 
directly to 8 otW = V x U (Rogachevskii & Kleeorin 2003). 

There is evidence that the strong dynamo action seen in the simulations 
is only possible due to the combined presence of open boundaries and shear. 
This however has so far only been checked explicitly for the a effect that 
is present when the forcing is helical; see Sect. 15.61 In the case of the solar 




Importance of magnetic helicity in dynamos 19 




Fig. 9. Saturation behavior of the ratio q = (B )/(B 2 ) for runs with different 
kinetic helicity of the flow. Solid line: zero helicity, dotted line: positive helicity 
(opposite to the sun) dashed line: negative helicity (as in the sun). 



surface such losses are actually observed to occur in the form of coronal mass 
ejections and in active regions. In the sun, coronal mass ejections are quite 
vigorous events that are known to shed large amounts of helical magnetic 
fields (Berger & Ruzmaikin 2000, DeVore 2000, Chae 2000, Low 2001). This 
kind of physics is not at all represented by adopting vacuum or pseudo- 
vacuum (vertical field) boundary conditions, as was done in Brandenburg & 
Sandin (2004). 



5 Connection with the a effect 
5.1 Preliminary considerations 

The a effect formalism provides so far the only workable mathematical frame- 
work for describing the large scale dynamo action seen in simulations of 
helically forced turbulence. (In this section we retain the /xo factor.) The 
governing equation for the mean magnetic field is 

fyti 

— = Vx(UxB + £- Wo J) , (29) 



where £ = u x b is the electromotive force resulting from the u x b non- 
linearity in the averaged Ohm's law. Without mean flow, {7 = 0, and an 
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electromotive force given by a homogeneous isotropic a effect and turbulent 
diffusion r]t, i.e. 

£ = aB - ri t fj,oJ, (30) 

we have 

-=aVxS+( I , + t)t )V 2 B, (31) 

which has solutions of the form B = £} e lk x + xt w itfi the dispersion relation 

A± = -r lT k 2 ± \ak\, (32) 
and three possible eigenfunctions (appropriate for the periodic box) 



B(x) = 





or , (33) 



where k m = ki = 1. Obviously, when the coefficients a and tjt = ?7+??t remain 
constant, and there is an exponentially growing solution (for \a\ > ryxfci), the 
solution must eventually grow beyond any bound. At the latest when the 
magnetic field reaches equipartition with the kinetic energy, a and r/t must 
begin to depend on the magnetic field itself. However, the present case is 

2 

sufficiently simple so that we can continue to assume that B , as well as a 
and r]t, are uniform in space and depend only on time. 

Comparison with simulations has enabled us to eliminate a large number 
of various quenching models where a = a(B). The only quenching model 
that seems reasonably well compatible with simulations of a 2 -like dynamo 
action in a periodic box without shear is 

a = = 2 , rjt = zzo (empirical), (34) 

1 + R m B /£2 q i + RmB /B c 2 q 

see Fig- EI However, this type of quenching is not fully compatible with mag- 
netic helicity conservation, as has been shown by Field & Blackman (2002). 
This will be discussed in the next section. 



5.2 Dynamical a. quenching 

The basic idea is that magnetic helicity conservation must be obeyed, but the 
presence of an a effect leads to magnetic helicity of the mean field which has 
to be balanced by magnetic helicity of the fluctuating field. This magnetic 
helicity of the fluctuating (small scale) field must be of opposite sign to that 
of the mean (large scale) field. 

We begin with the uncurled mean-field induction equation, written in the 
form 

— = S - Tjfi J, (35) 
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dot it with B, add the result to A ■ dB/dt, average over the periodic box, 
and obtain 

j t (A-B) = 2(e ■B)-2r,n (J-B). (36) 

To satisfy the helicity equation for the full field, (A ■ B) = (A ■ B) + (a ■ b), 
we must have 

^- t (a-b) = -2(S-B)-2 mo (j-b). (37) 

Note the minus sign in front of the 2(£ ■ B) term, indicating once again that 
the a effect produces magnetic helicity of opposite sign at the mean and 
fluctuating fields. The sum of the two equations yields Eq. JBJ. 

The significance of Eq. (|37() is that it contains the (j ■ b) term which 
contributes to the a effect, as was first shown by Pouquet et al. (1976). 
Specifically, they found (see also Blackman & Field 2002) 

a = «K + Om, with qr = — |t(u; • u), om = +\ T {j ■ b), (38) 

where r is the correlation time of the turbulence, u) = V x it is the vorticity, 
and (uj ■ u) is the kinematic helicity. 

Using (j ■ b) w fc| (a • b), see Eq. (|13|) . we can rewrite Eq. (|37|) in a form 
that can directly be used in mean field calculations: 

da 2 / a{B 2 ) — rj t (J ■ B) a- a K \ 

* = ~ 2vtkt { — w q — + ) ' (39) 

Here we have used % = ^Tu^ ms to eliminate r in favor of rjt and B^ q = 
/ioPou^ ms to eliminate u^ ms in favor of B% . 

So, a is no longer just an algebraic function of B, but it is related to B via 
a dynamical, explicitly time-dependent equation. In the context of dynamos 
in periodic domains, where magnetic helicity conservation is particularly im- 
portant, the time dependence of a can hardly be ignored, unless one wants 
to describe the final stationary state, which can be at the end of a very slow 
saturation phase. However, in order to make contact with earlier work, it is 
useful to consider the stationary limit of Eq. i.e. set da/dt. 



5.3 Steady limit and its limitations 

In the steady limit the term in brackets in Eq. (|39|l can be set to zero, so this 
equation reduces to 

R ^ a(B )-n t (j-B) + ft = (for da/dt = Q y (4Q) 

Solving this equation for a yields (Kleeorin & Ruzmaikin 1982, Gruzinov & 
Diamond 1994) 
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a K + r h R m (J-B)/Bl l 
a= (tor da/at — 0). (41) 

l + R m (B )/B 2 q 

And, sure enough, for the numerical experiments with an imposed large scale 
field over the scale of the box (Cattaneo & Hughes 1996), where B is spa- 
tially uniform and therefore J = 0, one recovers the 'catastrophic' quenching 
formula, 

a= ^ (for 7 = 0), (42) 

l + R m {B )/ J B| q 

which implies that a becomes quenched when (B )/B 2 q = R^ 1 ~ 10~ 8 for 
the sun, and for even smaller fields in the case of galaxies. 

On the other hand, if the mean field is not imposed, but maintained by 
dynamo action, B cannot be spatially uniform and then J is finite. In the 
case of a Beltrami field 1)33(1 . (J -B)/(B ) = fc m is some effective wavenumber 
of the large scale field [k m = e m fc m ; see Eq. Since i? m enters both the 

numerator and the denominator, a tends to ry t fc m , i.e. 

a -> r/ t k m (for 7^0 and 7 || B). (43) 

Compared with the kinematic estimate, qk ~ f}t^ii a is only quenched by 
the modified scale separation ratio. More importantly, a is quenched to a 
value that is just slightly above the critical value for the onset of dynamo 
action, a cr it = r]^k m . How is it then possible that the fit formula l|34|) for a 
and rjt produced reasonable agreement with the simulations? The reason is 
that in the simple case of an a 2 dynamo the solutions are degenerate in the 
sense that J and B are parallel to each other. Therefore, the term (J ■ B)B 

is the same as (B ) J, which means that in the mean EMF the term aB, 
where a is given by Eq. I|41l) , has a component that can be expressed as being 
parallel to J. In other words, the roles of turbulent diffusion (proportional to 
J) and a effect (proportional to B) cannot be disentangled. This is the force- 
free degeneracy of a 2 dynamos in a periodic box (BB02) . This degeneracy is 
also the reason why for a 2 dynamos the late saturation behavior can also 
be described by an algebraic (non-dynamical, but catastrophic) quenching 
formula proportional to 1/(1 + R m (B 2 )) for both a and rjt, as was done in 
B01. To see this, substitute the steady state quenching expression for a, from 
Eq. l|41Jl . into the expression for £. We find 

S f , UK + RnMJ-B)/B 2 q 
xB - (rj + rft) J = ; — _ 2 — - B - rjt J 



l + R m (B*)/B 2 q 
a K B 7ft J 



l + R m (B 2 )/B 2 l + R ni (B 2 )/B 2 



(44) 



which shows that in the force-free case the adiabatic approximation, together 
with constant (unquenched) turbulent magnetic diffusivity, becomes equal to 
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the pair of expressions where both a and rj t are catastrophically quenched. 
This force-free degeneracy is lifted in cases with shear or when the large scale 
field is no longer fully helical (e.g. in a nonperiodic domain, and in particular 
in the presence of open boundaries). 



5.4 The Keinigs relation and its relevance 

Applying Eq. (|3"?|) to the steady state using £ = aB — r/tfioJ (and retaining 
/xo factor), we get 

a = — nun^— 9— + %k m (for periodic domain), (45) 
( B ) 

where we have defined an effective wavenumber of the large scale field, fc m = 
(io(J-B)/(B ) ; see Eq. © ■ This relation applies only to a closed or periodic 
box, because otherwise there would be boundary terms. Moreover, if the 
mean field is defined as a volume average, i.e. B = (B) = Bo, then ji^J = 
V x Bq = 0, so k m = and one has simply 

a = —T]fj, - (for imposed field). (46) 
( B ) 

This equation is due to Keinigs (1983). For the more general case with k m 7^ 
this equation has been discussed in more detail by Brandenburg & Sokoloff 
(2002) and Brandenburg & Matthaeus (2004). 

Let us now discuss the significance of this relation relative to Eq. (|41|l . 
Both equations apply only in the strictly steady state, of course. Since we 
have assumed stationarity, we can replace (j ■ b) by —(J-B); see Eq. ® . 
Thus, Eq. l(4*5|> reduces to 

a = 77T& m (47) 

where i]t = T) + rjt is the total (turbulent and microscopic) magnetic diffu- 
sivity. This relation is just the condition for a marginally excited dynamo; 
see Eq. (j22J), so it does not produce any independent estimate for the value 
of a. In particular, it does not provide a means of independently testing 
Eq. H41fl . The two can however be used to calculate the mean field energy in 
the saturated state and we find (BB02) 

2 

(B ) a — 77T&m s 

— 5— = (periodic domain). (48) 

#cq Vtkm 

By replacing k m by an effective value fc m , this equation can be generalized to 
apply also to the case with shear (for details see BB02). 
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5.5 Blackman's multi-scale model: application to helical 
turbulence with imposed field 

The restriction to a two scale model may in some cases turn out to be insuffi- 
cient to capture the variety of scales involved in astrophysical bodies. This is 
already important in the kinematic stage when the small scale dynamo obeys 
the Kazantsev (1968) scaling with a k 3 ^ 2 spectrum that peaks at the resistive 
scale. As the dynamo saturates, the peak moves to the forcing scale. This lead 
Blackman (2003) to develop a four scale model where he includes, in addition 
to the wavenumbers of the mean field k m (= k\) and the wavenumber of the 
energy carrying scale of the velocity fluctuations fcf (= ki), also the viscous 
wavenumber k v (= ks) and the resistive wavenumber k n (= £4). The set of 
helicity equations for the four different scales is 



where £\ is the usual electromotive force based on kinetic helicity at the 
forcing scale, &2, with feedback proportional to H 2 , and £ 2 has no kinetic 
helicity input but only reacts to the automatically generated magnetic helicity 
H3 produced at the viscous scale k%. These equations are constructed such 
that 



which is consistent with the magnetic helicity equation © for the total field. 
An important outcome of this model is that in the limit of large R m the 
magnetic peak travels from k 3 to k 2 on a dynamical timescale, i.e. a timescale 
that is independent of R m . 

Brandenburg & Matthaeus (2004) have applied the general idea to the 
case of a model with an applied field. Here the new scale is the scale of the 
applied field, but since this scale is infinite, this field is fixed and not itself 
subject to an evolution equation. Nevertheless, the electromotive force from 
this field acts as a sink on the next smaller scale with wavenumber k\ , which 
is the largest wavenumber in the domain of the simulation. They thus arrive 
at the following set of evolution equations, 



(d t + 2 V k 2 ) H x 
(d t + 2 V k 2 ) H 2 
(d t + 2r)kl) H 3 



2(£ 1 -B 1 )+2{£ 2 -B 1 ), 
-2(£ 1 -B 1 ) + 2(£ 2 -B 2 ), 
-2(£ 2 -B 1 )-2(£ 2 -B 2 ), 



(49) 
(50) 
(51) 
(52) 





(53) 




(54) 

(55) 
(56) 



(6 t + 2 v k 2 ) H Y 



2(£ -B )+2{£ 1 -B 1 ), 




The square brackets around the first equation indicate that this equation 
is not explicitly included. From the second equation Q55JI one can see that 
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Fig. 10. Evolution of the total magnetic helicity, H — Hi + Hi, as a function of t 
for different values of Bo , as obtained from the three-dimensional simulation. Note 
the change of sign at Bo ~ B, w 0.07. [Adapted from Brandenburg & Matthaeus 



there is a competition between two opposing effects: the a effect operating 
on the imposed field Bo and the a effect operating on the B\ field on the 
scale of the box. When the imposed field exceeds a certain field strength, 
Bq > B*, the former will dominate, reversing the sign of the magnetic helicity 
at wavenumber k\ . This is actually seen in the simulations of helically forced 
turbulence with an imposed field Bq; see Fig. We return to this at the 
end of this section. 

The work of Brandenburg & Matthaeus (2004) was motivated by earlier 
work of Montgomery et al. (2002) and Milano et al. (2003) who showed that, 
if the imposed magnetic field is weak or absent, there is a strong nonlocal 
transfer of magnetic helicity and magnetic energy from the forcing scale to 
larger scales. This leads eventually to the accumulation of magnetic energy 
at the scale of the box (Meneguzzi et al. 1981 , Balsara & Pouquet 1999, B01). 
As the strength of the imposed field (wavenumber k = 0) is increased, the 
accumulation of magnetic energy at the scale of the box (k = 1) becomes 
more and more suppressed (Montgomery et al. 2002). 

In order to solve the model equations, we have to make some assumptions 
about the electromotive force operating at fco and k\ . The large scale magnetic 
helicity production from the a effect operating on the imposed field is £ q ■ 
Bq = cvi-Bq. On the other hand, S\ at wavenumber k\ is given by 



(2004).] 



Si = a t Bx - rjtHoJt. 



(57) 
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To calculate (£ x ■ Bi) in Eqs and we dot Eq. lj57|) with B x , volume 
average, and note that Ho(Ji ■ -Bi) = fci-Bi and {B\) — k\\H±\. The latter 
relation assumes that the field at wavenumber ki is fully helical, but that it 
can have either sign. Thus, we have 

{£i-B 1 )=a i k 1 \H 1 \-r h klHi. (58) 

The a effects on the two scales are proportional to the residual magnetic 
helicity of Pouquet et al. (1976); see Eq. I|38|) . In terms of H 1 and H 2 = Hf 
we write 

ai = a K + \rk\Hx, (59) 
a 2 = a K + lrkjH 2 , (60) 

for the a effect with feedback from Hi and H 2 , respectively. 

For finite values of Bo, the final value of Hi is particularly sensitive to the 
value of qk and turns out to be too large compared with the simulations. This 
disagreement with simulations is readily removed by taking into account that 
ax = — \t{u) ■ u) should itself be quenched when B becomes comparable to 
B cq . Thus, we write 

a K =a K0 /(l + Bl/B^), (61) 

which is a good approximation to more elaborate expressions (Riidiger & 
Kitchatinov 1993). We emphasize that this equation only applies to «k and 
is therefore distinct from Eqs. l|3"I)l. IpBty. or 

Under the assumption that the turbulence is fully helical, the critical 
value B* of the imposed field can be estimated by balancing the two terms 
on the right hand side of Eq. I|56l) and by approximating, a « r/ t kf and 
(j ■ b) w k { B^ q . This yields 

B 2 jB 2 eq ^ V /r, t = R-\ (62) 

where the last equality is to be understood as a definition of the magnetic 
Reynolds number, see also BB02. For Bo > B* the sign of the magnetic 
helicity is the same both at k = 1 and at k = k{ , while for Bo < B* the signs 
are opposite. 

In Fig. El w e show the result of a numerical integration of Eqs l|55|l 
and i|56|l . Both the three-dimensional simulation and the two-scale model 
show a similar value of Bq « 0.06. ..0.07, above which Hi changes sign. This 
confirms the validity of our estimate of the critical value B* obtained from 
Eq. ||fj5J. Secondly, the time evolution is slow when B < B» and faster 
when Bo > B». In the simulation, however, the field attains its final level 
for Bo > B* almost instantaneously, which is not the case in the model. The 
significance of this discrepancy remains unclear. Nevertheless, the level of 
agreement between the simulations and 3-scale model is surprising, suggest- 
ing that the approach can indeed be quite useful. 
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Fig. 11. Evolution of magnetic helicity as a function of t for different values of Bo, 
as obtained from the two-scale model. [Adapted from Brandenburg & Matthaeus 
(2004).] 



5.6 Alpha effect with open boundaries and shear 

In a recent paper, Brandenburg & Sandin (2004) have carried out a range 
of simulations for different values of the magnetic Reynolds number, R m = 
u r ms/(i]k{), for both open and closed boundary conditions using the geometry 
depicted on the right hand panel of Fig. In order to measure a, a uniform 
magnetic field, Bq = const, is imposed, and the magnetic field is now written 
as B = Bq + V x A. They have determined a by measuring the turbulent 
electromotive force, and hence a = (£) ■ Bq/Bq. Similar investigations have 
been done before both for forced turbulence (Cattaneo & Hughes 1996, B01) 
and for convective turbulence (Brandenburg et al. 1990, Ossendrijver et al. 
2001). 

As expected, a is negative when the helicity of the forcing is positive, and 
a changes sign when the helicity of the forcing changes sign. The magnitudes 
of a are however different in the two cases: |cv| is larger when the helicity of 
the forcing is negative. In the sun, this corresponds to the sign of helicity in 
the northern hemisphere in the upper parts of the convection zone. This is 
here the relevant case, because the differential rotation pattern of the present 
model also corresponds to the northern hemisphere. 

There is a striking difference between the cases with open and closed 
boundaries which becomes particularly clear when comparing the averaged 
values of a for different magnetic Reynolds numbers; see Fig. El With closed 
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Fig. 12. Dependence of |(«}|/urms on J? m for open and closed boundaries. The 
case with open boundaries and negative helicity is shown as a dashed line. Note 
that for R m « 30 the a effect is about 30 times smaller when the boundaries are 
closed. The dotted line gives the result with open boundaries but no shear. The 
vertical lines indicate the range obtained by calculating a using only the first and 
second half of the time interval. [Adapted from Brandenburg & Sandin (2004).] 

boundaries a tends to zero like , while with open boundaries a shows 
no such decline. There is also a clear difference between the cases with and 
without shear together with open boundaries in both cases. In the absence of 
shear (dashed line in Fig. ll2[) a declines with increasing R m , even though for 
small values of R m it is larger than with shear. The difference between open 
and closed boundaries will now be discussed in terms of a current helicity 
flux through the two open open boundaries of the domain. 

5.7 Current helicity flux 

It is suggestive to interpret the above results in terms of the dynamical a 
quenching model. However, Eq. (|39f) has to be generalized to take the diver- 
gence of the flux into account. In order to avoid problems with the gauge, 
it is advantageous to work directly with j ■ b instead of a ■ b. Using the evo- 
lution equation, db/dt = — V x e, for the fluctuating magnetic field, where 
e = E — E is the small scale electric field and E = r/fiQ J — £ the mean 
electric field, one can derive the equation 




= -2e~c 




(63) 



where 




2e X j + (V X e) X b/Mo, 



(64) 
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is the current helicity flux from the small scale field, and c = V x j the 
curl of the small scale current density, j = J — J. In the isotropic case, 
e • c « kfe-b, where k{ is the typical wavenumber of the fluctuations, here 
assumed to be the forcing wavenumber. Ignoring the effect of the mean flow 
on £ [as is usually done; but see Krause & Radler (1980) and the recent on the 
shear current effect by Rogachevskii & Kleeorin (2003, 2004); see Sect. 14.5] . 
we obtain 



-(« x B Q ) ■ b + rj/ioj b = £ B + r}fi j • b, (65) 



where we have used u x b = £ and Bo = B. Using standard expressions for 

3' 

strength, B eq = ^fiopu rms , we eliminate r via 



the turbulent magnetic diffusivity, r] t = hru 2 ms , and the equipartition field 



\t = noPoVt/B^. (66) 
This leads to an explicitly time dependent formula for a, 

This equation is similar to that of Kleeorin et al. (2000, 2002, 2003) who 
considered the flux of magnetic helicity instead of current helicity. 

Making use of the adiabatic approximation, i.e. putting the rhs of Eq. I|67ll 
to zero, one arrives at the algebraic steady state quenching formula (da/dt = 

°) 

a K + Rm {vmJ ■ B — |fc f 2 V • fio^W) /B 2 CCi 

a — — ^ • l""J 

1 + R m B /BZ q 

In the absence of a mean current, e.g. if the mean field is defined as an average 
over the whole box, then B = Bo = const, and J = 0, so Eq. (|68|l reduces to 

a K - \kfR m V ■ Mo^| S /< /Rn , 

a = — — -. (69) 

1 + R m B /B eq 

This expression applies to the present case, because we consider only the sta- 
tistically steady state and we also define the mean field as a volume average. 

For closed boundaries, (V-^f) = 0' anc ^ so Eq- l|6l)|l clearly reduces to a 
catastrophic quenching formula, i.e. a vanishes in the limit of large magnetic 
Reynolds numbers as 

^(closed) = ^ R -l ( for Rm ^ ( ?0) 

1 + R m B /B eq 

The i?" 1 dependence is confirmed by the simulations (compare with the 
dash-dotted line in Fig. I12|l . On the other hand, for open boundaries the 
limit i? m — > oo gives 
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a (o P c„) ^ _ (v . w ,^|S)/(2fc f 2 Bg) (for R m - oo), (71) 

which shows that losses of negative helicity, as observed in the northern hemi- 
sphere of the sun, would enhance a positive a effect (Kleeorin et al. 2000). 
In the simulations, the current helicity flux is found to be independent of 
the magnetic Reynolds number. This explains why the a effect no longer 
shows the catastrophic i?" 1 dependence (see Fig. lT2|) . In principle it is even 
conceivable that with qk — a current helicity flux can be generated, for 
example by shear, and that this flux divergence could drive a dynamo, as 
was suggested by Vishniac & Cho (2001). It is clear, however, that for finite 
values of i? m this would be a non-kinematic effect requiring the presence of 

1/2 

an already finite field (at least of the order of _B cq /i?m ). This is because 
of the 1 + R^B^/ B 2 ^ term in the denominator of Eq. I|69(l . At the moment 
we cannot say whether this is perhaps the effect leading to the nonhelically 
forced turbulent dynamo discussed in Sect. 14.51 or whether it is perhaps the 
S x J or shear-current effect that was also mentioned in that section. 



6 What about rj quenching? 

As we have seen above, in a closed domain the value of a in the saturated state 
cannot conclusively be determined without also determining at the same time 
the turbulent magnetic diffusivity. There are different ways of determining rjt . 
The values are not necessarily all in agreement with each other, because one 
does not know whether the mean field equation, where % enters, is correct 
and applicable. We report here a few different examples where rjt has been 
determined. 



6.1 Direct measurements in a working dynamo 

We first consider the case of a helical turbulent dynamo without shear (B01) 
and compare it with a simple mean-field a 2 dynamo. Assuming that a is 
uniform, we can use Eq. I|31[) and, assuming furthermore that a < (which 
is the case when the helicity of the forcing is positive, as in B01), the solution 
is 

B = (b x cos/ciz, by sinfciz, 0) T . (72) 
The time-dependent equations can then be written as 

= \a\b y -r lT k 2 b x , (73) 

^ = \a\b x - r, T k 2 b y . (74) 
In an isotropic, homogeneous a 2 dynamo, the eigenfunction obeys b x — b y . 
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Fig. 13. Result for rjr for different values of i? m . The lines represent the fits 
described in the text. In the plot of tjt the asterisks denote \a\ — A for the i? m ,f or c = 
120 run, which agrees reasonably well with r/T- [Adapted from B01.] 



We now assume that, at some particular time, we put b x = 0, for exam- 
ple. This means that b x (t) will first grow linearly in time at a rate that is 
proportional to a like b x w |a|fci& a . At the same time as b x grows, b y will 
first decrease at a rate that is proportional to 7]t- This allows an independent 
estimate of b x and b y by solving the matrix equation 

kxby -kfb x \ (\a\\ fdb x /dt\ 

hb x -kfbyj \r lT ) {dby/dtj- { ' 0) 

The result for a is found to be roughly consistent with that of Cattaneo & 
Hughes (1996), and the result for rjr is reproduced in Fig. ^3 and can be 
described by the fit formula 

* = 1 + ~ 9 \B\/B C(1 (?6) 

with g w 16. This expression needs to be compared with that obtained from 
other approaches. 

The fact that the results obtained for a by using this approach are con- 
sistent with that for a uniform field is quite surprising and unexpected. This 
agreement probably indicates that in this type of simulation a is independent 
of scale - at least in the scale range corresponding to wavenumbers k = k% 
(= 1) and k = 0. In general, this may not be true. Indeed, in the case of 
accretion discs some numerical evidence for scale dependence of a and rjt has 
been found (Brandenburg & Sokoloff 2002). 
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Table 1. Summary of the main properties of the three-dimensional simulations 
with shear. Here, rj/{c s k\) is the magnetic diffusivity in units of the sounds speed 
and the wavenumber of the domain, and ui cyc = 2ir/T cyc is the cycle frequency. In 
Run (iii) there is no clear cycle visible. [Adapted from BDS02.] 



Run 


(i) 


(") 


(iii) 




io- a 


5 x 10~ 4 


2 x 10" 


u/ri 


5 


10 


25 


flm = (M 2 ) 1/2 /(7?fcl) 


30 


80 


200 


Cn = (U 2 ) 1/2 /ivki) 


1000 


2000 


4000 




4 


6 


20 


20 


30 


60 


e a = no(J ■ B)/{B 2 ) 


0.11 


0.06 


0.014 


Ucyc /(Vk'i) 


8. ..9 


6. . . 12 


> 10? 



6.2 Measurements in an ex ft dynamo 

In the case of an afl dynamo the cycle frequency uj cyc depends directly on 
the nonlinearly suppressed value of ryx 

w cyc = VTk 2 (vt is quenched), (77) 

see BB02 (their Sect. 4.2). The estimates of BBS01 indicated that the dy- 
namo numbers based on shear, Cq = S/(r]Tkj), is between 40 and 80, whilst 
the total dynamo number (V — C a Cn) is between 10 and 20 (see BBS01), 
and hence C a = a/(r/Tki) ~ 0.25. Thus, shear dominates strongly over the a 
effect {Cn/Ca is between 150 and 300), which is typical for aJ7-type behav- 
ior (i.e. oscillations) rather than a 2 -type behavior which would start when 
Co/C a is below about 10 (e.g. Roberts & Stix 1972). 

The results shown in Table ^ suggest that the period in this oscilla- 
tory dynamo is controlled by the microscopic magnetic diffusivity because 
w cyc /(r]kl) is approximately independent of R m . Using Eq. H77fl . this means 
that r] T (quenched)/ 7? = O(10) for R m between 30 and 200. This result would 
favor a model where ?7t is still quenched in an i? m -dependent fashion. In 
the next section we show that the apparent i? m -dependent r/ t quenching can 
easily also be produced when the field possesses a helical component. 

Looking at the scaling of the cycle frequency with resistivity may be quite 
misleading in the present case, because the large scale magnetic field exceeds 
the kinetic energy by a large factor (20-30). This would always lead to the 
usual (non-catastrophic) quenching of a and r)t- Furthermore, such strong 
magnetic fields will affect the mean shear flow. Most important is perhaps 
the fact that in the simulation of BBS01 the shear flow varies sinusoidally in 
the cross stream direction, so the mean field depends on the two coordinate 
directions perpendicular to the streamwise direction. For this reason BB02 
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Table 2. Results from the simulations of BBS01 and BDS02, compared with those 
of 2-dimensional mean field models. Model results that are in fair agreement with 
the simulations are highlighted in bold face. Here, Q is the ratio of toroidal to 
poloidal rms field. 



Model 






C n 


9 


S 

■qk' 2 


(b 2 ) 


(B 2 ) 
B 2 


Q- 1 


e m 


^cyc 

~~S~ 


A 
S 


BBS01 


80 


1-2 






2000 


6 


30 


0.014 


0.06 


0.008 


0.015 


Rl 


20 


1.0 


100 





2000 


0.20 


15 


0.031 


0.065 


0.016 


0.044 


AG 2 


100 


0.5 


20 


3 


2000 


0.10 


22 


0.011 


0.024 


0.006 


0.021 


BDS02 


30 


1-2 






1000 


4 


20 


0.018 


0.11 


0.014 


0.006 


s3 


30 


0.35 


33 


1 


1000 


0.07 


6 


0.029 


0.061 


0.014 


0.016 


SI 


30 


0.35 


33 


3 


1000 


0.07 


19 


0.009 


0.019 


0.005 


0.016 



solved the mean field and dynamical quenching equations in a 2-dimensional 
model. It turned out to be important to allow for non-catastrophic quenching 
of rjt using Eq. i|76|) where the value of g has been varied between and 
3. The asymptotic l/B behavior (as opposed to 1/B 2 , for example) was 
motivated both by simulations (B01) and analytic results (Kitchatinov et al. 
1994, Rogachevskii & Kleeorin 2001). 

In order to see whether the models can be made to match the direct 
simulations, several input parameters were varied. It should be kept in mind, 
however, that not all input parameters are well known. This has to do with the 
uncertainty in the correspondence between the magnetic Reynolds number 
in the model (which measures r/to /v) an d the simulations [where it is defined 
as u rms / (r/k{)]. Likewise, the dynamo number C a = a/{rfrk\) is not well 
determined. Nevertheless, many of the output parameters are reasonably well 
reproduced; see Table 

6.3 Decay experiments 

Finally, we consider the decay of a magnetic field. This provides a fairly 
straightforward method of determining ?7t from the decay rate A of a sinu- 
soidal field with wavenumber kx, so A = r]^k\. The result reported by Yousef 
et al. (2003) suggests that 

v t *i h = (0.8 . . . 0.9) x u ims /k { (for B 2 « B c 2 q ). (78) 

Once the mean flow profile has decreased below a certain level (< 0.1w rms ), 
it cannot decay further and continues to fluctuate around 0.08u rms , corre- 
sponding to the level of the rms velocity of the (forced!) turbulence at k = k% 
(see the dashed line in Fig.[T3)l. 

The quenching of the magnetic diffusivity, r}t — rjt(B), can be obtained 
from one and the same run by simply determining the decay rate, Xb(B), 
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Fig. 14. Decay of large scale helical velocity and magnetic fields (dashed and 
solid lines, respectively). The graph of U(t) has been shifted so that both U(t) 
and B(t) share the same tangent (dash-dotted line), whose slope corresponds to 
ft = rjt — 0.86u rms /fcf. The decay of a nonhelical magnetic field is shown for 
comparison (dotted line). [Adapted from Yousef et al. (2003)] 

at different times, corresponding to different values of B = \B\. To describe 
departures from purely exponential decay one can adopt a Independent rj t 
expression of the form (|?FJ)l . It turns out that the value of g is not universal 
and depends on the field geometry. This is easily demonstrated by comparing 
the decay of helical and nonhelical initial fields; see Fig. ED 

In the next section we show that the slower decay of £?, and hence the 
implied stronger quenching of i] t , can also be described by a self-induced 
magnetic a effect which acts such as to decrease the decay rate. In the case 
of a helical initial field, we have J x B = 0, i.e. the large scale field is force-free 
and interacts only weakly with the turbulence. 

Thus, the indications here are that for non-helical fields, rjt is not catas- 
trophically quenched. A resistively slow decay rate occurs however when the 
magnetic field is helical, but this is not to be explained by a catastrophically 
quenched n t , but by the magnetic a effect, that tries to keep the magnetic 
field as large as possible, just as enforced by the magnetic helicity constraint. 
The phenomenon, described in this way, may be more easily described in 
terms of helicity conservation, because the system has magnetic helicity that 
can only decay slowly on a resistive time scale, hence lowering the apparent 
turbulent diffusivity down to the microscopic value r\. This will be explained 
in more detail in the next section. 
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Fig. 15. Dependence of the turbulent diffusion coefficient on the magnitude of the 
mean field. R m « 20. Left: The initial field is helical and corresponds to data points 
on the right hand side of the plot. The data are best fitted by g — 8 = 0.4_R m . Right: 
the same for the nonhelical case. The data are best fitted by a — I, independent of 
R m . [Adapted from Yousef et al. (2003)] 



6.4 Taylor relaxation or selective decay 

In the case of a helical field with B /B^ q <; i?" 1 the slow decay of B is related 
to the conservation of magnetic helicity. As already discussed by BB02, this 
behavior is related to the phenomenon of selective decay (e.g. Montgomery 
et al. 1978) and can be described by the dynamical quenching model. This 
model applies even to the case where the turbulence is nonhelical and where 
there is no a effect in the usual sense. However, the magnetic contribution to 
a is still non-vanishing, because it is driven by the helicity of the large scale 
field. 

To demonstrate this quantitatively, Yousef et al. (2003) have adopted the 
one-mode approximation (fc = fei) with B = Bexp(ifciz), the mean-field 
induction equation 

^ = ifer x £ - nkfB, (79) 
together with the dynamical a-quenching formula (|39|l . 



da „ , 2 / - Rc(5 • B 
— = -2 V kf a + i? m ^^ 



2r)kf \ a + -Rm^T^ L , (80) 

oq / 



where 

£ = aB - Tftifci x B (81) 

is the electromotive force, and i? m is defined as the ratio rjto/r], which is 
expected to be close to the value of R m . 
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Fig. 16. Dynamical quenching model with helical and nonhelical initial fields. 
The quenching parameters are 5 = (solid line) and 3 (dotted line) . The graph for 
the nonhelical cases has been shifted in t so that one sees that the decay rates are 
asymptotically equal at late times. The value of ifr used to normalize the abscissa 
is based on the unquenched value. [Adapted from Yousef et al. (2003)] 

In Fig. E| we show the evolution of B/B cq for helical and nonhelical 
initial conditions, B oc (1, i, 0) and B oc (1,0,0), respectively. In the case of 
a nonhelical field, the decay rate is not quenched at all, but in the helical 
case quenching sets in for B /B^ q ^ . In the helical case, the onset 

of quenching at B /B% q ss i?" 1 is well reproduced by the simulations. In 
the nonhelical case, however, some weaker form of quenching sets in when 
B /-Bg q pa 1 (see the right hand panel of Fig, IT^|l . We refer to this as standard 
quenching (e.g. Kitchatinov et al. 1994) which is known to be always present; 
see Eq. (|76ll . BB02 found that, for a range of different values of R m , g — 3 
resulted in a good description of the simulations of cyclic ai?-type dynamos 
(BDS02). 

Yousef et al. (2003) also showed that the turbulent magnetic Prandtl 
number is indeed independent of the microscopic magnetic Prandtl number. 
The resulting values of the flow Reynolds number, Re = u rms / (vkf), varied 
between 20 and 150, giving P m in the range between 0.1 and 1. Within plot 
accuracy the three values of As turn out to be identical in the interval where 
the decay is exponential. 

7 Conclusions 

In the present review we have tried to highlight some of the recent discoveries 
that have led to remarkable advances in the theory of mean field dynamos. 
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Of particular importance are the detailed confirmations of various aspects of 
mean field theory using helically forced turbulence simulations. The case of 
homogeneous turbulence with closed or periodic boundary conditions is now 
fairly well understood. In all other cases, however, the flux of current helicity 
becomes important. The closure theory of these fluxes is still a matter of on- 
going research (Kleeorin et al. 2000, 2002, 2003), Vishniac & Cho (2001), Sub- 
ramanian & Brandenburg (2004), and Brandenburg & Subramanian (2004). 
The helicity flux of Vishniac & Cho (2001) has been independently confirmed 
(Subramanian & Brandenburg 2004). A more detailed investigation of current 
helicity fluxes appears to be quite important when one tries to get qualitative 
and quantitative agreement between simulations and theory. 

The presence of current helicity fluxes is particularly important when 
there is also shear. This was already recognized by Vishniac & Cho (2001) 
who applied their calculations to the case of accretion discs where shear is 
particularly strong. In the near future it should be possible to investigate the 
emergence of current helicity flux in more detail. This would be particularly 
interesting in view of the many observations of coronal mass ejections that 
are known to be associated with significant losses of magnetic helicity and 
hence also of current helicity (DeVore 2000, Demoulin et al. 2002, Gibson et 
al. 2002). 

In order to be able to model coronal mass ejections it should be par- 
ticularly important to relax the restrictions imposed by the vertical field 
conditions employed in the simulations of Brandenburg & Sandin (2004). A 
plausible way of doing this would be to include a simplified version of a corona 
with enhanced temperature and hence decreased density, making this region 
a low-beta plasma. 

In the context of accretion discs the importance of adding a corona is 
well recognized (Miller & Stone 2000), although its influence on large scale 
dynamo action is still quite open. Regarding hydromagnetic turbulence in 
galaxies, most simulations to date do not address the question of dynamo 
action (Korpi et al. 1999, de Avillez & Mac Low 2002). This is simply be- 
cause here the turbulence is driven by supernova explosions which leads to 
strong shocks. These in turn require large numerical diffusion, so the effec- 
tive magnetic Reynolds number is probably fairly small and dynamo action 
may only be marginally possible. In nonhelically driven turbulence has been 
applied to the galactic medium to argue that it is dominated by small scale 
fields (Schekochihin et al. 2002), but the relative importance of small scale 
fields remains still an open question (Haugen et al. 2003). Galaxies are how- 
ever rotating and vertically stratified, so the flows should be helical, but in 
order to say anything about magnetic helicity evolution, much larger mag- 
netic Reynolds numbers are necessary. At the level of mean field theory the 
importance of magnetic helicity fluxes is well recognized. The explicitly time- 
dependent dynamical a quenching equation with magnetic helicity fluxes has 
been included in mean field simulations (Kleeorin et al. 2000, 2002, 2003), 
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but the form of the adopted fluxes is to be clarified in view of the differences 
with the results of Vishniac & Cho (2001) and Subramanian & Branden- 
burg (2004). Nevertheless, given that the form of the dynamical quenching 
equations is likely to be still incomplete, it remains to be demonstrated, us- 
ing simulations, that magnetic or current helicity fluxes do really allow the 
dynamo to saturate on a dynamical time scale. 
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